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Abstract. In this paper, we give several simple methods for drawing a whole rational 
surface (without base points) as several Bezier patches. The first two methods apply to 
surfaces specified by triangular control nets and partition the real projective plane RP 2 
into four and six triangles respectively. The third method applies to surfaces specified by 
rectangular control nets and partitions the torus RP 1 x RP 1 into four rectangular regions. 
In all cases, the new control nets are obtained by sign flipping and permutation of indices 
from the original control net. The proofs that these formulae are correct involve very little 
computations and instead exploit the geometry of the parameter space (MP 2 or RP 1 x MP 1 ). 
We illustrate our method on some classical examples. We also propose a new method for 
resolving base points using a simple "blowing up" technique involving the computation of 
"resolved" control nets. 



1 Introduction 



In this paper, we consider the problem of drawing a whole rational surface. For example, 
consider the sphere F specified by the fractions 

. 2u 2v . u 2 + v 2 — 1 

x{u,V) = 2 , 2 _,-, , y{u,V) = — 2——, Z{U,V) = — 

U 2 + V 2 + 1 + f 2 + 1 U 2 + V 2 + 1 

The problem is that no matter how large the interval [r, s] is, the trace F([r,s] x [r, s]) of 
F over [r, s] x [r, s] is not the trace of the entire surface. In this particular example, we 
could take advantage of symmetries, but in general, this may not be possible. We could 
use any of the bijections from ] — 1, 1[ to K. to reduce the parameter domain to the square 
[— 1, 1] x [—1, 1], but since these maps are at least quadratic, this could triple the total degree 
of the surface, leading to an impractical method. For example, using the map 



1 -t 2 

,u f 1 u (l-u 2 ) 2 (l-v 2 ) 
the traction — becomes 



u 2 + uv (1 — v 2 )u 2 + (1 — u 2 )uv 

Recomputing the control net after substitution would also be quite expensive. Indeed, one 
of the reasons why the problem is not trivial is that in most CAGD applications, the surface 
is given in terms of control points rather than parametrically (in terms of polynomials). 

Thus, the problem is to cope with the situation in which u or v become infinite. But what 
do we mean exactly by that? To deal with this situation rigorously, we can "go projective", 
that is, homogenize the polynomials. However, this can be done in two different ways. The 
first method is to homogenize with respect to the total degree, replacing u by u/t and v by 
v/t, getting 

2ut _ 2vt u 2 + v 2 - 1 2 

U 2 + V 2 + t 2 ' ^ U 2 + V 2 + t 2 ' U 2 + V 2 + t 2 

The parameter domain is now the real projective plane MP 2 . Points at infinity are the points 
of homogeneous coordinates (u, v,0) (i.e., when t = 0). Observe that all these point at 
infinity yield the north pole (0, 0, 1) on the sphere. 

The second method is to homogenize separately in u and v, replacing u by u/t 1 and v by 
v/t 2 , getting 

2uM! _ 2vt\t 2 _ u 2 t\ + v 2 t\ - t\t\ 

X ~^ 2 T^W+M' V ~ u 2 t 2 + v 2 t 2 + t 2 t 21 Z ~ uH\ + vH\ + t\t\ ' 

This time, the parameter domain is the product space MP 1 x MP 1 , where MP 1 is the real 
projective line. The domain MP 1 x MP 1 is homeomorphic to a torus, and it is not homeo- 
morphic to RP 2 . Observe that when ti — t 2 — 0, all the numerators and the denominator 
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vanish simultaneously. We have what is called a base point. This is annoying but not terribly 
surprising, since a sphere is not of the same topological type as a torus. It should be noted 
that there are also rational surfaces (such as the torus) that do not have base points when 
treated as surfaces with domain RP 1 x RP 1 but have base points when treated as having 
domain MP 2 (see Section 6), and vice versa. 

In summary, there are two ways to deal with infinite values of the parameters. We can 
homogenize with respect to the total degree m (replacing u by u/t and v by v/t). This leads 
to rational surfaces specified by triangular control nets, as we will see more precisely in the 
next section. The other method is to homogenize with respect to u and the maximum degree 
p in u (replacing u by u/ti) and with respect to v and the maximum degree q in v (replacing 
v by v/t?). This leads to rational surfaces specified by rectangular control nets, as we will 
see more precisely in the next section. 

The problem of drawing a rational surface reduces to the problem of partitioning the 
parameter domain into simple connected regions Ri such as triangles or rectangles, in such a 
way that there is some prespecified region R and some projectivities such that every other 
region is the image of the region R Q under one of the projectivities. Furthermore, if the 
patch associated with the region R is given by a control net Ao, we want the control net 
Mi associated with the region Ri to be computable very quickly from J\f . 

In the case of the real projective plane RP 2 , we can use the fact that RP 2 is obtained 
as the quotient of the sphere S 2 after identification of antipodal points. The real projective 
plane can be partitioned by projecting any polyhedron inscribed in the sphere S 2 on a plane. 
This way of dividing the real projective plane into regions is discussed quite extensively in 
Hilbert and Cohn-Vossen [14] (see Chapter III). As noted by Hilbert, it is better to use 
polyhedra with central symmetry, so that the projective plane is covered only once since 
vertices come in pairs of antipodal points. In particular, we can use the four Platonic solids 
other than the tetrahedron, but if we want rectangular or triangular regions, only the cube, 
the octahedron, and the icosahedron can be used. Indeed, projection of the dodecahedron 
yields pentagonal regions (see Hilbert and Cohn-Vossen [14], page 147-150). 

If we project the cube onto one of its faces from its center, we get three rectangular 
regions (see Section 3). It is easy to find the projectivities that map the central region onto 
the other two. Since we are dealing with the projective plane, it is better to use triangular 
control nets to avoid base points, and it is necessary to split the central rectangle into two 
triangles. Thus, the trace of the rational surface is the union of six patches over various 
triangles. It is shown in Section 3 how the control nets of the other four patches are easily 
(and cheaply) obtained from the control nets of the two central triangles. 

If we project the octahedron onto one of its faces from its center, we get four triangular 
regions (see Section 4). This time, it is a little harder to write down the projectivities that 
map the central triangle rst to the other three triangles R, S, T. However it is not necessary 
to find explicit formulae for these projectivities, and using a geometric argument, we can find 
very simple formulae to compute the control nets associated with the other three triangles 
from the control net associated with the central triangle, as shown in Section 4. 
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Projecting the icosahedron onto one of its faces from its center yields ten triangular 
regions, but we haven't found formulae to compute control nets of the other regions from the 
central triangle. We leave the discovery of such formulae as an open and possibly challenging 
problem. 

Let us now consider the problem of partitioning PJP 1 x PJP 1 into simple regions. Since the 
projective line PJP 1 is topologically a circle, a very simple method is to inscribe a rectangle 
(or a square) in the circle and then project it. One way to do so leads to a partition of MP 1 
into [—1, 1] and RP 1 — ] — 1, 1[. The corresponding projectivity is t \— > j. Other projections 
lead to a partition of MP 1 into [r, s] and MP 1 — ]r, s[ for any affine frame (r, s). In all cases, the 
torus is split into four rectangular regions, and there are very simple formulae for computing 
the control nets of the other three rectangular nets from the control net associated with the 
patch over [—1, 1] x [—1, 1] (or more generally, [r±, Si] x [r 2 , s 2 ]), as explained in Section 5. 

It should be stressed that it is not necessary to compute explicitly the various projec- 
tivities, and that in each case, a simple geometric argument yields the desired formulae for 
the new control nets. Other methods for drawing rational surfaces were also investigated by 
Bajaj and Royappa [2, 3] and DeRose [5] and will be discussed in Section 4 and Section 5. 

There is a problem with our methods when all the numerators and the denominator 
vanish simultaneously. In this case, we have what is called a base point. In Section 6, we 
give a new method for resolving base points (in the case of a rational surface specified by a 
triangular control net), using a simple "blowing-up" technique based on an idea of Warren 
[18]. What is new is that we give formulae for computing "resolved" control nets. 

It turns out that to give rigorous proofs of our formulae, it is necessary to view rational 
surfaces as surfaces defined in a suitable projective space in terms of multiprojective maps. 
We will summarize how to do this in Section 2. The proofs that our formulae are correct 
involve very little computations and instead exploit the geometry of the parameter space 
(MP 2 or MP 1 x MP 1 ). For the sake of brevity, we do not review how polynomial surfaces 
are defined in terms of control points. The deep reason why polynomial triangular surface 
patches can be effectively handled in terms of control points is that multivariate polynomials 
arise from multiaffine symmetric maps (see Ramshaw [17], Farin [7, 6], Hoschek and Lasser 
[15], or Gallier [11]). 

2 Rational Surfaces and Control Points 

Denoting the affine plane M 2 as V, a rational surface F : V — > W 1 of degree m is specified by 
some fractions 



where Fi(u, v), . . . , F n+ i{u, v) are polynomials of total degree < m. In order to handle ra- 
tional surfaces in terms of control points, it turns out that it is necessary to view rational 
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surfaces as surfaces in some projective space. Roughly, this means that we have to homog- 
enize the polynomials Fi(u, v), . . . , F n+ i(u, v). However, the polar forms of homogeneous 
polynomials are multilinear, and thus we must deal with multilinear maps rather than mul- 
tiaffine maps. Fortunately, there is a construction to embed an affine space into a vector 
space, in such a way that multiaffine maps extend uniquely to multilinear maps. This con- 
struction is described in Berger [4] and is at the heart of the presentation of rational surfaces 
in Fiorot and Jeannin [8, 9]. However, Fiorot and Jeannin do not use polar forms. We have 
adapted their approach in the framework of polar forms in Gallier [12]. In this paper, we 
simply review the facts needed to understand the proof of our theorems. Given a vector 
space E, we denote the projective space induced by E as F(E) (see Berger [4] or Gallier 
[12]). Given an affine space E with associated vector space E , a vector space E can be con- 
structed, such that E is embedded as an affine hyperplane in E via an affine map j : E — > E, 
and E as a hyperplane. 
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Figure 1: Embedding an affine space into a vector space 

Both hyperplanes are defined by some linear form ou: E — > R. The previous diagram 
illustrates the embedding of the affine space E into the vector space E: Roughly, the vector 
space E has the property that for any vector space F and any affine map / : E — > F , there 
is a unique linear map f : E —> F extending / : E —> F . As a consequence, given two affine 
spaces E and F, every affine map / : E — > F extends uniquely to a linear map / : E — > F. 

A pair (a, A) where a G E and A ^ is called a weighted point. Vectors in E of the form 
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(a, 1) are identified with points a in E. It is easily shown that for every a G E, we have 
E = ~E © Ra. 

We have the following important result whose proof can be found in Gallier [12], or 
inferred from Ramshaw [17]. 

Lemma 2.1 Given any two affine spaces E and F and a multiaffine map f: E m — > F, there 
is a unique multilinear map /: (E) m — > F extending f as in the diagram below: 

E m f * F 

jx-xj j 

{E) m —^ F. 

Given an affine space E, the projective space F(E) induced by E is denoted as E, and it 
is called the projective completion of E. Observe that R = 1RP 1 and V = MP 2 . 

The upshot of the above considerations is that a rational surface can be defined in terms 
of multilinear maps. Let £ be some ambiant affine space in which our surfaces live, in most 
cases M 3 . If we first homogenize the polynomials Fi(u,v) with respect to the total degree m 
(replacing ubj u/z and v by v/z), we can view a rational surface as a map 

F: V -> £ 

(where V = M. 2 is the affine plane) induced by some symmetric multilinear map 

/: (V) m ^£ 

such that 

F({u,v,z]) = nf)((^v,z),...,(u,v,z)), 

" v ' 

m 

for all homogeneous coordinates (u,v,z) e M 3 . We call such surfaces rational total degree 
surfaces, or triangular rational surfaces. Furthermore, for any affine frame Arst, the trian- 
gular control net Af = (&i,j,k)(i,j,k)eA m (in £) w.r.t. Arst defining the triangular surface F is 
given by the formulae 

Oi,j,k = /Qv— ^_-j_5> jj;), where i + j + k = m. 

i j k 

On the other hand, if we first homogenize the polynomials Fi(u,v) with respect to u 
and the maximum degree p in u (replacing u by u/ti) and second with respect to v and the 
maximum degree q in v (replacing v by v/t 2 ), we can view a rational surface as a map 

F:RxR^£ 
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induced by some multilinear map 

/: (§) p x (R) q 

which is symmetric in its first p arguments and in its last q arguments, and such that 
F([u, tx], [v, t 2 ]) = P(/)((u, h), . . . , (u, h), (v, i 2 ), . . . , (v, t 2 )), 

V v ' N * ' 

p 9 

for all homogeneous coordinates (u,ti), (v,t 2 ) G IR 2 . We call such surfaces rational surfaces 
of bidegree (p, q), or rectangular rational surfaces. Furthermore, given any two afline frames 
(ri, si) and (r2, S2) for the affine line R, the rectangular control net M = {9i,j)o<i<p,o<j<q ( m 
£) w.r.t. (ri, Si) and (r 2 , s 2 ) defining the rectangular surface F is given by the formulae 

hi = /C r i» • • • , si, . . . , s 1? r 2 , . . . , r 2 , s 2 , . . . , s 2 ^ 



p—i 



1-3 



A 6ase point of a rational surface F specified by a multilinear map / is any point a G 
such that f(a, . . . , a) = , or any point (w, t> ) G M 2 such that f(u, . . . , u, v, . . . , v) — . 



3 Splitting Triangular Rational Surfaces Into Six Tri- 
angular Patches 

As we explained in Section 1, if we project a cube onto one of its faces from its center, we 
obtain a partition of the projective plane MP 2 into three rectangular regions, in such a way 
that there exist simple projectivities f and ip between the square [—1, 1] x [—1, 1] and the 
other two regions. 



V 




u 



Figure 2: Dividing the projective plane into 3 rectangular regions 
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The projectivity is induced by the linear isomorphism of IR 3 given by 

(u, v, w) I— > (v, w, u). 

Choosing the line at infinity w — 0, the restriction of this map to the afline plane V (cor- 
responding to w — 1) is the map (u,v, 1) i— > (v, l,u). This is the map from V to V such 
that (u,v) i — (f/w, l/-u) if u ^ 0, and (0, v) \— > (t>, 1,0) when w = 0. The projectivity ^ is 
induced by the linear isomorphism of M 3 given by 

(u, v, w) I— > (ty, -u, f ). 

Choosing the line at infinity w — 0, the restriction of this map to the afline plane P (corre- 
sponding to w — 1) is the map (u, i>, 1) (l,-u,f). This is the map from V to V such that 
(u, v) t— > (1/f , -u/t>) if f ^ 0, and (m, 0) i— > (1, -u, 0) when f = 0. 

Actually, it turns out that the method of this section holds for any region defined by a 
nondegenerate quadrilateral (a,b,c,d), i.e. when (a,b,c,d) is a projective frame. However, 
the details are a bit messy, and for simplicity, we restrict out attention to a rectangular region 
[ri, si] x [r 2 , s 2 ]. Since we are dealing with triangular surfaces, it will be necessary to split the 
rectangle [n, si] x [r 2 , s 2 ] into two triangles, and thus, we will obtain the trace of a rational 
surface as the union of 6 patches over various triangles in the rectangle [ri,Si] x [r 2 ,s 2 ]. 
Letting a, b, c, d be the vertices of the rectangle [ri, si] x [r 2 , s 2 ] defined such that 

a = (si,s 2 ), b=(n,s 2 ), c=(ri,r 2 ), d=( Sl ,r 2 ), 

as shown below 



v 



b 


a 


















c 


d 



Figure 3: Some afline frames associated with the rectangle (a, b, c, d) 
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we will consider the following afflne frames 



Abca = ((r u s 2 ), (ri,r 2 ), (si,s 2 )), 
Aetac = ((si,r 2 ), (si,s 2 ), (?"i,r 2 )), 
A&ad = ((ri,s 2 ), (si,s 2 ), (si,r 2 )). 

In particular, a rectangular surface patch defined over the rectangle [n, si] x [r 2 , s 2 ] will 
be treated as the union of two triangular surface patches defined over the triangles Abca and 
Adac. It is somewhat unfortunate that a control net over the third frame Abad needs to be 
computed, but that is what the proof of lemma 3.2 shows. In any case, such a control net 
can be computed very cheaply from a control net over Abca (or Adac). 

There is simple geometric explanation of the partitioning method in terms of the usual 
model of the real projective plane V = MP 2 in M 3 . Recall that in this model, the real projec- 
tive plane MP 2 consists of the points in the plane z — 1 corresponding to the lines through 
the origin not in the plane z — 0, and of the points at infinity corresponding to the lines 
through the origin in the plane z — 0. We view the vertices of the rectangle (a, b, c, d) defined 
above as points in the plane z — 1, in which case their coordinates are (s±, s 2 , 1), (ri, s 2 , 1), 
(ri, r 2 , 1), and (si, r 2 , 1). Then, we have the parallelepiped (a, b, c, d, —a, —6, — c, — d). There 
is a unique projectivity P(y?) such that 

P(v?)(a) = a, F(<p)(b) = c, P(</?)(c)=d, P(^)(d) = 6. 

For instance, it is induced by the unique linear map </? such that 

ip(a) = a, <f(b) = — c, <f(c) = —d. 

Since d = — b + a + c, we get 

</?(d) = -</?(&) + ^{ a ) + v( c ) = c + a — d — ( r i> s 2, 1) = 

The linear map (p transforms the top face (a, b, c, d) of the parallelepiped to the back face 
(a, — c, —d,b). When a line L through the origin and passing through a point of the face 
(a, — c, — d, b) varies, the intersection of L with the plane z — 1 varies in ty?([ri, Si] x [r 2 , s 2 ]). 
We can define a rhombus (a, e, /, g, —a, — e, — /, — g>) inscribed in the sphere of center O = 
(0, 0, 0) and of radius R = a/s 2 + s 2 + 1 passing through a, as follows: the points e, /, g 
are on the upper half-sphere and they are determined by the intersection of the lines (O, b), 
(O, c) and (O, d) with the sphere. Then, it is obvious that under the central projection of 
center O onto the plane z — 1, the top face (a, e, /, g) of the rhombus projects onto the face 
(a, b, c, <i) of the parallelepiped, and that the projection of the rhombus onto the plane z — 1 
yields the desired partitioning of MP 2 . 

Similarly, there is a unique projectivity P(^) such that 

P(V>)(a) = a, P(V0( & ) = p W0(c) = 6, P(V0(d) = c. 
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It is induced by the unique linear map if) such that 

if>(a) = a, if)(b) = d, if>(c) = —b. 

Since d = — b + a + c, we get 

if)(d) = —if)(b) + if) (a) + if)(c) = —d + a - b = (— n, — r 2 , -1) = — c. 

The linear map ^ transforms the top face (a, 6, c, cf) of the parallelepiped to the right face 
(a,d, —b, — c). When a line L through the origin and passing through a point of the face 
(a, d, —b, — c) varies, the intersection of L with the plane z = 1 varies in -0 ( [r"i , Si] x [r 2 , S2]). 
Again, it is obvious that under the central projection of center O onto the plane z = 1, 
the top face (a, e, /, g) of the rhombus projects onto the face (a, b, c, d) of the parallelepiped, 
and that the projection of the rhombus onto the plane z = 1 yields the desired partition- 
ing: of MP 2 . Fi gure 4 shows the parallelepiped (a, b, c, d, —a, —b, — c, —d) and the rhombus 
(a,e,f,g,-a,-e,-f,-g). 




Figure 4: Parallelepiped and rhombus associated with (a, b, c, d) 

We will now use the maps (p and if) to show how the trace of a rational surface F 
can be obtained as the union of the traces of three rational surfaces over the rectangle 
[r 1; si] x [r 2 , S2]. 1 The first of these surfaces is F itself, and the two other rational surfaces 

^^While reading Appell's Treatise of Rational Mechanics, we stumbled on the fact that the change of 
variable (u, v) ^ (l/v,u/v) was used by Appell in his solution to a problem of Bertrand (see [1], Tome I, 
Part III, Chapter XI, page 422-423). Appell explains that he found this "homographic transformation" in 
1889. The problem of Bertrand is to find all central force laws depending only on the position of a moving 
particle, so that the trajectory of the particle is a conic for every choice of initial conditions. 
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and F$ are easily obtained from F. However, depending on the multilinear map / defining 
F, the surface F (and thus, F v and F^) may have base points, that is, we may have 

f{ j.u,v,z),...,(u,v,z)) = ~0 

m 

for some (u, v, z) ^ (0, 0, 0). We will show how to deal with this situation later on. 

In order to render the trace of F, we will use the fact that it is the union of the six traces 
F(Abca), F(Adac), F(ip(Abca)), F(ip(Adac)), F(ip(Abca)), and F(ip(Adac)). Furthermore, 
the last four traces are also obtained as traces of F v and F$ over some appropriate choice of 
affine frames among Abca, Adac, and Abad. 

We now show how F v and F^ are defined, and how their control points can be computed 
very simply from the control points of F (computed with respect to the affine frames Abca, 
Adac, and Abad). We will assume that the homogenization V of the affine plane V is 
identified with the direct sum IR 2 © M.O, where O = (0,0). Then, every element of V is of 
the form (u, v, z) G M 3 . 

Definition 3.1 Given an affine space £ of dimension > 3, for every rational surface F: V — > 
£ of degree m specified by some symmetric multilinear map /: (V) m — > £, the symmetric 
multilinear maps f v : (V) m — > £ and \ (V) m — > £ are defined such that 

U((ui, «i, wi),..., (u m , v m , w m )) = f((p(ux, v x , w x ),..., ip(u m , v m , w m )), 

U((U1, Vx, Wx),..., (u m , V m , W m )) = f{lp{llx, Vx, Wx),..-, 1p(u m , V m , W m )). 

Let Ftp-. V — > £ be the rational surface specified by f^\ (V) m — > £, and let F^: V — > £ be 
the rational surface specified by fy: (V) m — * £. 

Observe that the base points of F v , if any, have coordinates (u, v, w) ^ (0, 0, 0) such that 

f((p(u, v, w), . . . , ip(u, v, w)) = , 

and that the base points of F^, if any, have coordinates (u,v,w) ^ (0,0,0) such that 

f(ip(u, v, w), . . . , ip(u, v, w)) = . 

Lemma 3.2 Given an affine space £ of dimension > 3, for every rational surface F:V^>£ 
of degree m specified by some symmetric multilinear map f: (V) m — > £, if f 9 and are 
the symmetric multilinear maps of definition 3.1, except for base points, F , F^ and F^ have 
the same trace. The trace of F v over Abca is the trace of F over ip(Abca), the trace of F^ 
over Adac is the trace of F over ip(Adac), the trace of F^ over Abca is the trace of F over 
tp(Abca), and the trace of F^ over Adac is the trace of F over ip (Adac) . Furthermore, if 
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the control nets (in S) of the surface F w.r.t. the affine frames Abca, Adac, and Abad, are 
respectively 

ol = (ai,j,k)(i,j,k)eA m , 

ft = (Pi,j,k)(i,j,k)eA m , 
7 = {li,j,k)(i,j,k)€A rn , 

the control nets 9 1 and 9 2 (in S) of the surface F v w.r.t. the affine frames Abca and Adac, 
and the control nets p 1 and p 2 (in S) of the surface w.r.t. the affine frame Abca and 
Adac, are given by the equations 

0}, j ,k = (-l) i+i Pi,k,i, 
^i,* = (-l) fc 7i,i,fc, 
P\,i,h = {- l ) 3 lj,k,i, 

Proof. We have 

U((ui, Vi, Wi),..., (u m , v m , w m )) = f(tp(u!, Vi, Wi),..., ip{u m , v m , w m )), 

and thus 

P(/<^)(K,^l,Wl], • • • , [u m ,V m ,W m ]) = ¥(f)([<p(u 1 ,Vi,W 1 )], . . ., [(p(Um,V m ,W m )]). 

In view of the properties of </?, it is clear that F and F v have the same trace (except for base 
points), and that the trace of F v over Abca is the trace of F over ip(Abca), and the trace 
of F,p over Adac is the trace of F over ip(Adac). A similar argument applies to F and F^. 
The formulae for computing the control points of F v w.r.t. the triangle Abca are obtained 
by computing 

f v Qv^A c^^c, a^^a) . 

i j k 

Since 



U((U1, V U Wi),..., (Um, V m , W m )) = f((p(u!, V U IWi), . . . , (p(u m , V m , W m )), 

(fib) = — c, ip(c) = —d, and if (a) = a, we have 

f v (b, ...,b,c,...,c,a,..., a) = f(-c, -c, -d, . . . , -d, a, . . . , a), 



that is 



f v (b, . . . ,b,c^^,a^^a) = (-l) t+3 f (c^^, d^^d,, a^^a), 

j k i j k 
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and since the control points f3ij,k are computed w.r.t. the triangle Adac, we get 

(- \y j h.,,,. 

The formulae for computing the control points of w.r.t. the triangle Adac are obtained 
by computing 

f v ( d, ...,d , a,..., a , c,...,c) . 

i j k 

Since 

/ v ((ui, V U Wi),..., (u m , v m , w m )) = f((p(u 1: Vi, «7i), . . . , ^(u ro , u m , w m )), 
</?(d) = 6, <f(c) = —d, and </?(a) = a, we have 

/^(d, . . . , d, a, . . . , a, c, . . . , c) = f(b, . . . , b, a, . . . , a, — d, . . . , — d), 

j k i j k 



that is 



/Jd, . . . , d, a, . . . , a, c, . . . , c) = (-l) k f(b, . . . , b, a, . . . , a, d, . . . , d), 



and since the control points 7i,j,fe are computed w.r.t. the triangle Abad, we get 

^i, fc = (-l) fc 7i,;,*- 

The formulae for computing the control points of w.r.t. the triangle Abca are obtained 
by computing 

/y, <v^c, a, . . . , a) . 

i j k 

Since 

/^((Ul, Ui, Wi), • • • , (?Z m , U m , W m )) = f(lp(u u V U Wx), . . . , ^(Mm, Um, W m )), 

■0(6) = d, ^(c) = —6, and ^(a) = a, we have 

/^(6, . . . , 6, c, . . . , c, a, . . . , a) = /(d, . . . , d, -b, . . . , -b, a, . . . , a), 



that is 



f^(b, . . . , b, c, . . . , c, a, . . . , a) = . . . , d, b, . . . , b, a, . . . , a), 



and since the control points jij t k are computed w.r.t. the triangle Abad, we get 

Pl, j: k = {-l) j l3,k,i- 
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Finally, the formulae for computing the control points of w.r.t. the triangle Adac are 
obtained by computing 

/y, (d, . , d, cv^a , c^^c) . 

i j k 

Since 

/^,((ui, Ui, wi), • • • , (w m , u m , w m )) = «i, • • • , ^(u m , u m , iu m )), 

= — c, ^(c) = —6, and ^(a) = a, we have 

/^(o^^,a^^ = /( r c, . , -c , a, . „ , a, -b, . . . , -& ), 

i j fc i j k 

that is 

i j k i j k 

and since the control points ol%^^ are computed w.r.t. the triangle Abca, we get 

= (-!)*** <W 

□ 

The above calculations show that </? and ^ can be defined as above provided that d = 
—b + a + c, or equivalently b + d = a + c, which means that (a, b, c, d) is a parallelogram. 
Actually, lemma 3.2 also holds in the more general situation where (a, b, c, d) is a projective 
frame, i.e. a quadrilateral whose vertices are in general position. However, the definition of 
the linear maps tp and ip is a little more messy. As before, we identify a, b, c, d with points in 
the plane z — 1, and we let a = (a 1; a 2 , 1), 6 = (&i, 62, 1), c — (c±, c 2 , 1), and d = (di, d 2 , 1). 
To find a linear map ip inducing the unique projectivity P(v?) such that 

%)(«) = «, PH(6) = c, P(^)(c)=d, P(p)(d) = &, 

we let d = \a + fib + vc and b = X'a + fx'c + v'd, where A + /i + v — 1 and A' + // + v' — 1, 
and is the unique linear map such that 

ip(\a) = X'a, <p(fJ>b) = fjfc, <p(vc) = v'd. 

Then, <p(d) = b, as desired. The linear map ip can be defined in a similar way. The proof 
still goes through since the maps involved are multilinear, and thus not disturbed by scalar 
multiples. 

Lemma 3.2 shows that in order to render a rational surface, provided that it does not 
have base points, we just need to compute the control nets a, /3, 7 for the surface F w.r.t. 
the affine frames Abca, Adac, and Abad, since then, the control nets 6 l and 6 2 (in £) of the 
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surface w.r.t. the affine frames Abca and Adac, and the control nets p 1 and p 2 (in £) of 
the surface w.r.t. the affine frame Abca and Adac, are obtained at trivial cost. 

Remark: It should be noted that the surface patches associated with the control nets a, 
[3, 6 1 , 6 2 , p 1 , and p 2 , may overlap in more than boundaries. In fact, there are examples 
where a and (3 determine the entire surface, and other examples in which 6 1 , 6 2 , p 1 , and p 2 , 
determine the entire surface. 

It is fairly easy to implement this method in Mathematica. The interested reader will 
find such an implementation in Gallier [12]. In the interest of brevity, we content ourselves 
with some examples. 

Example 1. The algorithm is illustrated by the following example of an ellipsoid defined by 
the fractions 

. , 2ci« , 2c 2 v , c 3 (u 2 + v 2 -l) 

X ( U i V ) = —, 2~TT' y\ U ' V ) = 2 i 2 j 1 ' Z{U,V) = — ■ 

U 2 + V 2 + 1 U 2 + V 2 + 1 U 2 + V 2 + 1 

It is easily verified that this representation of the ellipsoid is derived from the stereo- 
graphic projection from the north pole onto the plane z — 0. The coordinates of a point on 
the sphere are the coordinates of the image of a point (u, v) the xOy plane, under the inverse 
of stereographic projection. We leave as an exercise to show that the following triangular 
control net for C\ = 4, c 2 = 3, c 3 = 2, is obtained: 

net = {{0, 0, -2, 1}, {0, 3, -2, 1} , {0, 3, 0, 2}, 
{A, 0, -2, 1}, {4, 3, -2, 1}, {4, 0, 0, 2}} 

The following picture shows the result of iterating the subdivision algorithm 3 times on 
the nets netl and net2: 
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4 

Figure 5: Patches 1, 2, of an ellipsoid 
Iterating the subdivision algorithm 3 times on the nets thetal and theta2 yields: 




Figure 6: Patches 3, 4, of an ellipsoid 
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Iterating the subdivision algorithm 3 times on the nets rhol and rho2 yields: 




Figure 7: Patches 5, 6, of an ellipsoid 
The result of putting all these patches together is the entire ellipsoid: 




4 



Figure 8: An entire ellipsoid 
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Of course, we could have taken advantage of symmetries, and our point is to illustrate 
the algorithm. 

Example 2. The Steiner roman surface is the surface of implicit equation 

x 2 y 2 + y 2 z 2 + x 2 z 2 = 2xyz. 

It is easily verified that the following parameterization works: 

2v . . 2u . . 2uv 

x{u, V) = - ■ — 5——, y(u,v) = ^— — 5——, z{u,v) 



M 2 + f 2 + l' ' u 2 + v 2 + l' u 2 + v 2 + 1 

It can be shown that this surface is contained inside the tetrahedron defined by the planes 

—x + y + z — 1, x — y + z — 1, x + y — z — 1, —x — y — z — 1, 

with —1 < x,y,z < 1. The surface touches these four planes along ellipses, and at the 
middle of the six edges of the tetrahedron, it has sharp edges. Furthermore, the surface is self- 
intersecting along the axes, and is has four closed chambers. A more extensive discussion can 
be found in Hilbert and Cohn-Vossen [14] , in particular, its relationship to the heptahedron. 
A triangular control net is easily obtained: 

steinl = {{0, 0, 0, 1}, {1, 0, 0, 1}, {1, 0, 0, 2}, 
{0, 1, 0, 1}, {1, 1, 1, 1}, {0, 1, 0, 2}}; 

We can display the entire surface using the method described in this section. Indeed, 
all six patches are needed to obtain the entire surface. One view of the surface obtained by 
subdividing 3 times is shown below (see Figure 9). Patches 1 and 2 are colored blue, patches 
3 and 4 are colored red, and patches 5 and 6 are colored green. A closer look reveals that 
the three colored patches are identical under appropriate rigid motions, and fit perfectly. 




Figure 9: The Steiner roman surface 
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Figure 10: A cut of the Steiner roman surface 



Another revealing view (see Figure 10) is obtained by cutting off a top portion of the 
surface. This way, it is clear that the surface has chambers. 

4 Splitting Triangular Rational Surfaces Into Four Tri- 
angular Patches 

As explained in Section 1, we obtain a partition of the real projective plane RP 2 into four 
triangles if we project an octahedron onto one of its faces from its center. We sketch such a 
method, leaving the simple details to the reader. 
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Figure 11: Splitting MP 2 into four triangles 

We let r, s, t be the vertices of the central triangle. The four triangles defined by the 
lines (r, s), (s,t), and (r,t) are denoted as rst, R, S, and T, where R,S,T contain points 
at infinity. It is easy to find three projectivities ipi\ MP 2 — > MP 2 , i = 1,2,3, such that 
ipi(rst) = R, (p2(rst) = S, and (p 3 (rst) = T. Then, we get some rational surfaces Fi = Fotpi, 
i — 1,2, 3. Indeed, if we use the model of MP 2 in R 3 where the points r, s, t are considered 
as being in the plane z — 1, it is immediately verified that the linear maps 

(r, s, t) I— > (— r, s, t), 
(r,s,t) ^ (r, -s,t), 
(r,s,t) ^ (r,s,-t), 

induce ipi,(p2i<P3- Furthermore, if the control net (in S) of the triangular surface F w.r.t. 
the affine frame Arst is a = (cxij,k)(i,j,k)eA m i ^ can be shown that the control nets 9 l ,9 2 , 
and 9 3 of the surfaces Fi, F 2 , F 3 w.r.t. Arts are given by the formulae 



1 


(-l)X 




2 


(-1)^ 




3 


(-1)% 





Provided that there are no base points, the traces of F,F 1 ,F 2 ,F 3 over Arst cover the 
entire trace of F (over MP 2 ). The upshot is that in order to draw a whole rational surface 
given by a triangular net a over Arst, we simply have to draw the four patches specified by 
a, O 1 ^ 2 , and 6* 3 , over Arst. 
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For example, we can apply the above method to the Steiner roman surface specified by the 
triangular net given in Example 2. It turns out that the patch F 3 is quite distorted. Applying 
the method to a net over a bigger triangle helps reduce the distorsion. In particular, we can 
send r and s to infinity, in which case the method ends up being equivalent to a method due 
to Bajaj and Royappa [2, 3]. Their method is based on the observation that the four maps 

/ <7\U a 2 v \ 

(u,v) ^ , , 

\1 — u — v 1 — u — V J 

where Oi G {—1, 1} for i — 1,2, map the triangle ((1, 0), (0, 1), (0, 0)) bijectively onto the four 
quadrants of the plane respectively. However, they do not consider the problem of computing 
the control nets of the surfaces 

\1 — u — v 1 — u — V J 

Another method for drawing triangular rational surfaces was also investigated by DeRose 
[5] who credits Patterson [16] for the original idea behind the method. Basically, the method 
consists in using the homogeneous Bernstein polynomials u l v^w k , where % + j + k = m, 
and to view a triangular rational surface as a rational map from the real projective plane. 
Then, by using any 3D model of the projective plane, it is possible to draw whole rational 
surface in one piece. For example, DeRose suggests to use an octahedron. However, the 
problem of finding efficient ways of computing control points is not addressed. 



5 Splitting Rectangular Rational Surfaces Into Four 
Rectangular Patches 

In this section, we show that every rectangular rational surface can be obtained as the union 
of four rectangular patches, and that the control nets for these patches can be computed 
very easily from the original control net. The idea is simple: we partition RP 1 x RP 1 into 
the four regions associated with the partitioning of MP 1 into [—1,1] and MP 1 — [—1,1]. Let 
(p be the projectivity of RP defined such that 

<p(u,t) = (t,u). 

We also define the following rectangular surfaces. 

Definition 5.1 Given an afnne space £ of dimension > 3, for every rectangular rational 
surface F: R x R — > £ of bidegree (p, q) specified by some (p, g)-symmetric multilinear map 
/: (R) p x (R) q — > £ , define the three (p, g)-symmetric multilinear maps fa: (M.) p x (M.) q — > £ , 
i — 1,2, 3, such that 



/l((ui, tl), . . . , (U p , t p ), (Ui, Si), . . . , (V q , S q )) 

/ 2 ((«i, h), . . . , (up, t p ), (vi, si), ... , {y q , s q )) 
fs((ui, h),..., (Up, t p ), (Vi, Si), ... , (v q , s q )) 



= f(<P(Ul, h),..., (p(Up, tp), (Vi, Si), ... , (Vq, S q )), 
= *l), • • • , (Up, tp), (f(Vl, Si), ... , ip(Vq, S q )), 

= f(<p(ui, h),..., tp(u p , tp), (pfa, Si), ... , <p(v q , Sq)). 
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The following lemma shows that provided that there are no base points, a rectangular 
rational surface is the union of four rectangular patches, and that given a rectangular net a 
w.r.t. (—1, 1) x (—1, 1), the other three nets can be obtained very easily from a. 

Lemma 5.2 Given an affine space £ of dimension > 3, for every rectangular rational surface 
F: R x R — > £ of bidegree (p, q) specified by some (p, q) -symmetric multilinear map f : (R) p x 
(R) 9 — > £, if fx, f 2 , fs are the {p, q) -symmetric multilinear maps of definition 5.1, except for 
the base points (if any), the trace 1, 1] x [—1, 1]) is the trace of F over ip([—l, 1]) x 

[— 1, 1], the trace F 2 ([— 1, 1] x [—1, 1]) is the trace of F over [—1, 1] x tp([— 1, 1]), and the trace 
F 3 ([— 1, 1] x [—1, 1]) is the trace of F over ip([—l, 1]) x (p([— 1, 1]). Furthermore, if the control 
net (in £) of the rectangular surface F w.r.t. (—1, 1) x (—1, 1) is 

a — ( a i,j)o<i<p,0<j<q, 

the control nets 9 1 , 9 2 , and 9 3 (in £) of the rectangular surfaces F 1: F 2: F 3 w.r.t. (—1, 1) x 
(—1,1) is are given by the equations 

(>:.., ( Q p ''>,.., 

91 = (-1)« 



a. 



a i,j- 



The proof is quite simple and left as an exercise. Actually, the same result applies to 
surfaces specified by a rectangular net over [r\, si] x [r 2 , s 2 ] for any affine frames (n, si) and 
(r 2 , s 2 ), since we can use the projectivity 

(s + r)t-2rs 
^ U 2t-(s + r) 

that maps [r,s] onto RP 1 — }r, s[. The upshot is that in order to draw a whole rational 
surface specified by a rectangular net a w.r.t. (ri, Si) x (r 2 , s 2 ), we simply have to compute 
the nets 9 1 , 9 2 , 9 3 , which is very cheap, and draw the corresponding rectangular patches. For 
example, a torus can be defined by the following rectangular net of bidegree (2, 2) w.r.t. 
(-1,1) x (-1,1): 

tornet4 = {{0 , -(a + b), 0, 4}, {0, 0, 4c, 0}, {0, (-a + b) , 0, 4}, 
{4(a + b), 0, 0, 0}, {0, 0, 0, 0}, {4(a - b) , 0, 0, 0}, 
{0, a + b, 0, 4}, {0, 0, 4c, 0}, {0, a - b, 0, 4}} 

For a = 2, b = 1, c = 1, we get 

tornet4 = {{0, -3, 0, 4}, {0, 0, 4, 0}, {0, -1, 0, 4}, 
{12, 0, 0, 0}, {0, 0, 0, 0}, {4, 0, 0, 0}, 
{0, 3, 0, 4}, {0, 0, 4, 0}, {0, 1, 0, 4}} 
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The result of subdividing the patches associated with F, Fi, F2 and F3 is shown below. 




Figure 12: A torus 



On the other hand, the method applied to a rectangular net of bidegree (2, 2) for an 
ellipsoid yields base points. For example, it can be shown that a control net of bidegree 
(2, 2) w.r.t. (—1, 1) x (—1, 1) for an ellipsoid is given by: 

recelnet3 = {{-8/3, -2, 2/3, 3}, {-8, 0, -2, 1}, {-8/3, 2, 2/3, 3}, 
{0, -6, -2, 1}, {0, 0, 6, -1}, {0, 6, -2, 1}, 
{8/3, -2, 2/3, 3}, {8, 0, -2, 1}, {8/3, 2, 2/3, 3}} 

Unfortunately, the patch corresponding to F 3 has a base point. The same thing happens 
for the Steiner roman surface. This is not surprising since neither the sphere nor the real 
projective plane are of the same topological type as the torus. It can be shown that a control 
net of bidegree (2, 2) w.r.t. (—1, 1) x (—1, 1) for the Steiner roman surface is given by: 

sqstein3 = {{-2/3, -2/3, 2/3, 3}, {0, -2, 0, 1}, {2/3, -2/3, -2/3, 3}, 
{-2, 0, 0, 1}, {0, 0, 0, -1}, {2, 0, 0, 1}, 
{-2/3, 2/3, -2/3, 3}, {0, 2, 0, 1}, {2/3, 2/3, 2/3, 3}}; 

Again, the the patch corresponding to F 3 has base points. 

Another method for drawing rectangular rational surfaces was investigated by DeRose 
[5] who credits Patterson [16] for the original idea behind the method. Basically, the method 
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consists in using the homogeneous Bernstein polynomials (?) (?) u l t\ % vH\ 1 , and to view a 
rectangular rational surface rational map from MP 1 x MP 1 . Then, by using any 2D 
model of the projective line, it is possible to draw a whole rational surface in one piece. 

In general, it is not easy to remove base points. This involves a technique from algebraic 
geometry known as "blowing-up" (see Fulton [10] or Harris [13]). In the next section, we 
will present a method for resolving base points in the case of triangular rational surfaces. 
However, we have not worked out the resolution of base points in the case of rectangular 
rational surfaces. We leave this problem as an interesting challenge to the reader. 

6 Resolving Base Points 

We now consider the case in which and (as defined in Section 3) have base points. An 
example for which this happens is the torus. 

Example 3. An elliptic torus can be defined parametrically as follows: 

x = (a — b sin <p) cos 9, 
y = (a — b sin ip) sin 9, 
z = c cos ip. 

As usual, we obtain a rational parameterization by expressing cost and sint in terms of 
tan(t/2), and we get the fractions 

(1 -u 2 ){a{l + v 2 ) -2bv) 
(l + u 2 )(l + v 2 ) ' 
2u(a(l + v 2 ) - 2bv) 
{l + u 2 ){l+v 2 ) ' 
c(l - v 2 ) 
1 + v 2 

Thus, the torus as a rational surface F is defined by 

x(u,v) = (1 -u 2 ){a{l + v 2 ) -2bv), 
y(u, v) = 2u(a(l + v 2 ) - 2bv), 
z(u,v) = c(l + u 2 )(l -v 2 ), 

w(u,v) = (1 + M 2 )(l + V 2 ). 

Rendering F over [—1,1] x [—1,1] yields one fourth of the torus, specifically, the front 
half of the upper half. Performing the change of variables 

(u,v) ^ ( -, - ), 
\u u J 



y = 

z = 
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the rational surface F v is defined by 

x(u,v) = (u 2 - v 2 )(a(l + u 2 ) -2bu), 
y(u, v) = 2uv(a(l + u 2 ) — 2bu), 
z(u,v) = c(u 2 - l)(u 2 + v 2 ), 
w(u,v) = (u 2 + v 2 )(u 2 + 1). 

Unfortunately, x(0, 0) = y(0, 0) = z(0, 0) = w(0, 0) = 0, and (0, 0) is a base point of F v . 
Performing the change of variables 



1 u 

t>' v 



(u,v) 

the rational surface F^ is defined by 

x(u, v) = (v 2 — l)(a(u 2 + v 2 ) — 2buv), 
y(u, v) = 2v(a(u 2 + v 2 ) — 2buv), 
z(u, v) = c{v 2 - u 2 ) (v 2 + 1), 
w(u,v) = (u 2 + v 2 )(v 2 + 1). 

Unfortunately, we also have x(0, 0) = y(0, 0) = z(0, 0) = w(0, 0) = 0, and (0, 0) is a base 
point of F^. 

If we try to render the rational surfaces F v and F^ over [—1, 1] x [—1, 1], we discover that 
some regions of these surfaces are not drawn properly. In these regions, there are holes and 
many lines segments shooting in all directions! The problem is that (0, 0) is a discontinuity 
point for both surfaces, and that the limit reached when u and v approach depends very 
much on the ratio v/u. One way to understand what happens is to let v = ku, simplify the 
fractions, and see what is the limit when u approaches 0. For F v , after calculations, we find 
that the limit when u approaches is 

a(l - k 2 ) 2ak 

i + k 2 ' TTP' 

which corresponds to the circle of radius a in the plane z = —c. For F^, after calculations, 
we find that the limit when u approaches is 

2bk n c(l - k 2 
-Q+ : , I9 , 0, 



1 + ' l + k 2 

which corresponds to an ellipse in the plane y — 0, centered at the point (—a, 0,0). It is 
indeed in the neighborhood of these two curves on the torus that F^ and F^ are not drawn 
properly. 
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We now propose a method to resolve the singularities caused by base points. The method 
is inspired by a technique in algebraic geometry known as "blowing-up" (see Fulton [10] or 
Harris [13]). What is new is that we give formulae for computing "resolved" control nets. 

In most cases, base points occur during a subdivision step in which a triangular net with 
a corner of zeros appears. Using a change of base triangle if necessary, it can be assumed 
without loss of generality that the corner of zeros has t as one of its vertices. If we display 
control nets (in £) with F(r) at the top corner, F(s) as the rightmost lower corner, and F(t) 
as the leftmost lower corner, a control net 9 = (Oij t k)(i,j,k)eA m of degree m has the following 
shape: 

x 

x x 

x x ... x x 
x x ... x x 
x x ... x x 

if ... if x x ... x x 
lflf...lflfxx...x x 

V v ' V v ' 

m m+l-rn 

It is assumed that all entries designated as x are nonzero. The more general case can be 
treated, but it is computationally too expensive to be practical. 

Given an affine frame Arst in the plane, recall that a rational surface F of degree m 
defined by the control net 9 = (Oi,j,k)(i,j,k)eA m is the projection onto £ of the polynomial 
surface G in £ defined by 9. Also, we have 

G(u,v) — J2 0i >* fc i^ uV ( 1-u ~ u )*' 

i+j+k=m 

for all u,v G KL It will be convenient to assume that if Qi,j,k G £^ is a weighted point, then 
its weight is denoted as Wij^, an d if 9i,j,k is a control vector, then we assign it the weight 
w i,j,k = 0. If we define w(u,v) as 



W(U,V)= W i,j,k^-WjU*V 3 (l-U-v) k , 



i+j+k=m 

whenever w(u,v) ^ 0, we have 

m! u l v J '(l — u — v) k 



F(u,v)= J2 °^ 



k 



i\j\k\ w(u,v) 

i+j+k=m J v ' ' 
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for all u, v G M. 

The "blowing-up" method used here relies on the following observation based on an idea 
of Warren [18]. Given the polynomial surface G in £ (and w), we define the polynomial 
surface G b and w b as follows: 

G b (a,P) = G(a(l-P),aP), 
Wb(ce,P) = w(a(l — ft), af3). 

Since a(l — 0) + a/9 = a, we get 

i+j+k=m ^ 

and 

= ^ w . J>fc ^ a *i(i _£)*. 

i+j+k=m 

Now, if ^i,j,fe = for i + j < to (with i + j + k = to), we note that both G b (a,P) and 
w b (a, P) are divisible by a m . If we define the polynomial surface G (and w), such that 

and 

w(a,p) = — , 

then we have _ 

G fc (a,/3) _ G(a,/3) 

w b (a,P) w{a,PY 
for all a ^ 0. Furthermore when a = 0, we have 

g(o,/?)= ^ m! ,. p j (i-py, 

ilv.lm — to! 

and 

™(o,/3)= v • T7-7— — — V] p j (i-py. 

ilv.lm — to! 

i+j=m J v ' 

Thus, for all /3 for which G(0, /3) and w(0,P) are not simultaneously null, 

g(M) 

w(0,/3) 
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is defined, and the polynomial surface G defines the rational surface F such that 

G(a,(3) 



F(a,0) = 



w(a, (3) 



Thus, what happens is that the triangular patch F over Arst is really a four-sided patch, 
the point F(t) being "blown up" into the rational curve of degree n whose control points are 

m—rn)i+j=m- 

If this rational curve^has no base points, then the rational surface patch F defined by the 
polynomial surface G has no base point, and it extends the surface patch F over Arst. If 
it has base points, they are common zeros of some polynomials in j3, and by simplifying by 
common factors and using continuity, we could eliminate these base points. For simplicity, 
we will assume that the boundary curve has no base points. 

Viewing G as a bipolynomial surface, note that G has bidegree (m — m, m). Also observe 
that the function 

{a, 13) h-> (a(l -P),a0) 
maps the unit square with vertices 

(0,0), (0,1), (1,1), (1,0) 

onto the triangle Arst = ((0, 1), (1,0), (0,0)), in such a way that the edge ((0,0), (0, 1)) is 
mapped onto t, the vertex (1, 1) is mapped onto s, and the vertex (1,0) is mapped onto r. 
Furthermore, if u = a(l — (3) and v = af3, we get 



a = u + v 



and 



P 



v 



U + V 



and thus, the map is invertible except on the line u + v = 0. Thus, we can think of the 
inverse map as "blowing up" the afnne frame Arst into the unit square. Specifically, the 
point t is "blown up" into the edge ((0, 0), (0, 1)). 



(0,1) 



(0,0) 



(1,0) 




Figure 13: Blowing up a triangle into a square 
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The only remaining problem is that the above method yields a rational square patch 
not given by a control net, and that we often need to map the unit square to an arbitrary 
base triangle. The second problem is easily solved. Assume that the afflne frame Arst has 
coordinates ( y (r 1 ,r 2 ), (s 1 , s 2 ), (ti,t 2 )). It is easily seen that the map defined such that 

u — (si - ri)aP + (r 1 - ti)(3 + t 1: 
v = (s 2 - r 2 )af3 + (r 2 - t 2 )(3 + t 2 , 

maps the unit square to the triangle Arst, in such a way that the edge ((0, 0), (0, 1)) is 
mapped onto t, the vertex (1, 1) is mapped onto s, and the vertex (1,0) is mapped onto r. 
Some simple calculations show that 

(s2 - r 2 )(u - ti) - (si - n)(v - t 2 ) 

a — 



P 



(ri - h)(s 2 - r 2 ) - (r 2 -t 2 )(si - n)' 
(ri - h)(v - t 2 ) - (r 2 - t 2 )(u - 



(s 2 - r 2 )(u - t x ) - (si - ri)(u - t 2 ) ' 
and thus, the map is only invertible outside the line of equation 

(s 2 - r 2 )(u - ij) - (si - r x )(v - t 2 ) = 0, 

the parallel to the vector (si — r\,s 2 — r 2 ) through t. 

Now, if g: (V) m — > £ is the polar form associated with G, we can compute the polar 
form g b : (V) m x (V) m — > £ associated with the bipolynomial surface G b as follows: 

g b (u X , . . . , U m , V 1: . . . , V m ) = — ^ 9{Ml - U ff (i)), • • • , («m(l - Uff(m)), «mU<r(m))), 

o"Go m 

where 5 m denotes the group of permutations on {1, ... , to}. The above formula corresponds 
to the case of the simple mapping u = a(l — (3), v = a/3, and it is obvious how to adapt it 
to the more general map 

u — (si — ri)aP + (ri - ti)@ + t 1: 
v = (s 2 - r 2 )a(3 + (r 2 - t 2 )/3 + t 2 , 

Now, over the afline basis (0,1), the square control net = (^□i,j)o<i,j<m associated 
with g b is defined such that 

0Oij = 9b(0,- ..,0,1,..., 1,0,. ..,0, 



m— i i m— j j 

However, if = for % + j < to, with % + j + k = to, then #Qj = for i < m, and thus 
we obtain the rectangular net 9Q of degree (to — to, to) associated with g, given by 



= (0Qj)m< 
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which corresponds to the rational surface defined by G. 

Thus, we know how to compute a rectangular net for the blown-up version G of G. A 
triangular net of degree 2m — to can easily be obtained. Indeed, there is a simple way for 
converting the polar form g: (V) p x (V) q — > £ of a bipolynomial surface of degree (p, q) into 
a symmetric multilinear polar form g& : (V) p+q — > using the following formula: letting 
to = p + g, we have 

9A{{ui,Vi),...,{u m ,v m )) = -4r f(n^n^)' 

Vp/ 7UJ={l,...,m} i6/ jGJ 

/nj=0 

|J|=P,|J|=9 

where 

II v i) = 9( u h, ■ ■ ■ ' Mi p' • • • ' 

with I = {i u . . .,i p }, and J = {ji, . . . ,j q }. 

Note that it is also possible to convert the polar form f\ '■ {V) m £ oi & surface of 
degree m into a symmetric (m, to) -multilinear polar form g: (V) m x (V) m — > ^, using the 
following formula: 

g(u u ...,u m ,v 1 ,...,v m ) = — ^2 • • • ; (V^m))), 

o"G<S m 

where S m denotes the group of permutations on {1, . . . , to}. 

Thus, we have a method for blowing up a control net 6 of degree to with a corner of zeros 
of size to, into a triangular net 6 of degree 2m — to, by first blowing up the triangular net 6 
into a rectangular net and then converting into a triangular net 0. 

Again, it is fairly easy to implement the above method in Mathematica (see Gallier [12]). 
In the interest of brevity, we content ouselves with some examples. 

Going back to Example 3 of this section, a torus, it turns out that in subdividing the nets 
thetal, theta2, rhol, and rho2, degenerate nets with a corner of zeros are encountered. 
In fact, these corners have two rows of zeros, and thus, the blowing up method yields nets 
of degree 6. For example, the net corresponding to thetal is resolved to a triangular net, 
which after 3 iterations of subdivision, yields the following picture: 
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Figure 14: A blow up of patch 3 of a torus 

The net corresponding to theta2 is resolved to a triangular net, which after 3 iterations 
of subdivision, yields the following picture: 




Figure 15: A blow up of patch 4 of a torus 
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Displaying these two pictures together, we get a shape reminicent of a horse-shoe crab! 




Figure 16: A blow up of patches 4,3 of a torus 
Similarly, blowing up the nets rhol and rho2 yields: 




Figure 17: A blow up of patches 5,6 of a torus 
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Together with the two patches associated with the square [—1, 1] x [—1, 1], we get the 
entire torus. Of course, we could have taken advantage of symmetries, and our point is to 
illustrate the algorithm. 

7 Conclusion 

We have presented several methods for drawing whole rational surfaces, defined paramet- 
rically or in terms of weighted control points. These methods rely on simple regular sub- 
divisions of the real projective plane or of the torus, and on versions of the de Casteljau 
algorithm. The main novelty is that the new control nets are obtained very cheaply from 
the original control net by sign flipping and permutation of indices. One of the advantages 
of our method is that it is incremental. Indeed, the algorithm produces an approximation 
of the surface as a sequence of control nets. Thus, if we wish to get better accuracy, we can 
subdivide each control net in the list. We can also achieve a zooming effect by selectively 
subdividing some subsequences of control nets. Bajaj and Royappa [2, 3] and DeRose [5] 
have also investigated methods for drawing whole rational surfaces. However, none of these 
papers address the problem of computing control nets. A weakness of our method is that it 
only applies to rational surfaces. On the other hand, although restricted to rational surfaces, 
our method is efficient, at least when there are no base points. We have also proposed a new 
method for resolving base points, by computing some refined control nets. Unfortunately, 
the present version of the method is exponential, and not practical as soon as the degree 
becomes greater than 4. Part of the problem is that our method first computes a rectan- 
gular control net which is then converted to a triangular net, and this conversion process is 
exponential. It would be interesting to compute directly a triangular net. 

Acknowledgement: We wish to thank Doug de Carlo for some very helpful comments. 
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